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Abstract. The X-ray afterglows of GRB001025A and GRB 010220 were detected by XMM-Newton with an average 0.2- 
10.0 keV flux of 4.4 and 3.3 x 10~ 14 ergcnr 2 s -1 respectively; the afterglow of GRB 00 1025 A is observed to decay. Afterglows 
at other wavelengths were not detected for either burst. A set of broadened soft X-ray emission lines are detected in the afterglow 
of GRB 001025 A, at 5.0 a significance above a Galactic-absorbed power-law continuum. The spectra of both afterglows are 
significantly better fit by a variable abundance thermal plasma model than by an absorbed power-law and are consistent with 
the observations of GRB 01 1211, indicating that thermal emission from light elements may be common in the early X-ray 
afterglows of GRBs. 
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1. Introduction 

Much of the recent progress in the understanding of gamma-ray 
bursts (GRBs) has come from bursts detected with good spatial 
accuracy with BeppoSAX and it is particularly at X-ray wave- 
lengths that GRB afterglows are detected (Piro 2001), about 
half producing no detectable optical afterglow emission (Fynbo 
et al. 2001). It is also only at X-ray wavelengths that emission 
lines are detected in afterglows, allowing firm estimates to be 
made of the cosmological redshifts and the outflow velocities 
of the afterglow material (Piro et al. 2000, Reeves et al. 2002). 
XMM-Newton (Jansen et al. 2001), with its large effective area, 
is particularly suited to this work. Previous detections of emis- 
sion lines in GRB afterglows with BeppoSAX and Chandra 
have concentrated on emission from highly-ionised iron (Piro 
et al. 1998, 2000, Antonelli et al. 2000, Yoshida et al. 2001); 
however recent observations with XMM-Newton have revealed 
several emission lines at lower energies (Reeves et al. 2002). 

Most plausible mechanisms for the production of a GRB 
involve a newly-formed black hole surrounded by a short-lived 
accretion disk regardless of the progenitor (Paczynski 1998, 
MacFadyen & Woosley 1999, Ruffert & Janka 1999, Meszaros 
2001). Recent evidence suggests that the progenitors of long- 
duration GRBs are massive stars (Paczynski 1998, Bloom et al. 
1999, MacFadyen & Woosley 1999, Reeves et al. 2002). 

Various models have been proposed to account for the 
emission spectra (in particular the claim of a high equivalent 
width Fe emission line (Antonelli et al. 2000)) and lightcurves 
of the afterglow. For instance, the nearby reprocessor model 
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(Rees & Meszaros 2000, Ballantyne & Ramirez-Ruiz 2001, 
Kallman et al. 2001) involving reflection of synchrotron emis- 
sion from the walls of a cone tunnelled out of a massive star, 
yielding large equivalent width Fe emission lines; or a 'supra- 
nova' model (Vietri & Stella 1998) which invokes a time de- 
lay (> 30 days) between an initial supernova (SN) explosion 
and the GRB, giving a spectrum briefly dominated either by 
the recombination of Fe in a photoionised plasma (Vietri & 
Stella 1998) or reflection of synchrotron emission off the walls 
of a wide funnel excavated in the SN remnant (Vietri et al. 
2001). Recently, Reeves et al. (2002) have suggested that the 
early X-ray afterglow spectrum of GRB 01 1211 is dominated 
by thermal emission from a metal-enriched, but notably Fe- 
poor collisionally-ionised plasma ejected in a recent SN explo- 
sion and heated by the GRB. 

In Sect. 2 we report on observations of two GRB afterglows 
with XMM-Newton, presenting the spectra in Sect. 3. In Sect. 4 
these results are discussed and their implications for other X- 
ray observations of afterglows examined. Our conclusions are 
in Sect. 5. Unless otherwise stated, all errors quoted are 90% 
confidence limits for one parameter of interest. A cosmology 
where Ho = 75 km s _I Mpc~' and qo = 0.5 is assumed through- 
out. 

2. Observations and data reduction 

2.1. GRB 001 025 A 

Smith et al. (2000) and Hurley et al. (2000) report detection 
of a GRB by the RXTE all-sky monitor, NEAR and Ulysses 
beginning at 03: 10:05 UT on 25 October 2000. It had a duration 
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of ~ 5 s, a fluence of 3.2 x 10~ 6 ergcirT 2 and a peak flux of 
2.0 x 1CT 6 erg crrT 2 s' 1 in the 25-100 keV band as detected by 
Ulysses. No variable source optical counterpart was detected 
down to an R-band magnitude limit of ~ 24.5 (Fynbo et al. 
2000). 

XMM-Newton began observing the error-box of 
GRB001025A 45 hours after the burst for two EPIC-pn 
exposures of 8ks and 16ks and single 37 ks exposures of the 
MOS cameras in full frame mode. The thin filters were used in 
each case. 

The data were processed and reduced with the SAS, version 
5.2, datasets from both EPIC-MOS cameras were co-added and 
both EPIC-pn exposures were also co-added and the resulting 
two datasets fit simultaneously. Source extraction regions were 
35" in radius and an off-source background extraction region of 
70" radius was chosen. Both single and double pattern (as well 
as triple and quadruple pattern for the MOS), good (FLAG=0) 
events were used with the ready-made response matrices 1 pro- 
vided by the XMM-Newton SOC. The spectra were binned with 
a minimum of 20 counts per bin. The final X-ray source posi- 
tions were determined after cross-correlation with the USNO 
A2.0 optical catalogue based on the SAS task eposcorr (see 
Tedds & Watson 2002). 

Two sources were detected in the IPN error-box (Altieri 
et al. 2000a,b); the brighter, with coordinates (J2000) RA. 
08h36m35.86s, Dec. -13°04'12.28", and a 68% error radius 
of 0.6", was seen to decay slowly (at 99.8% confidence), a 
power-law decay (F oc with index j3 = 3.0 ± 1.9 fitting 
the lightcurve well. We identify this source as the afterglow 
of GRB001025A. The Galactic hydrogen absorbing column 
in this direction is 6.1 x 10 20 crrr 2 (Dickey & Lockman 1990, 
using the FTOOL nh). Its mean 0.2-10.0keV flux was 4.4 x 
1 0" 14 erg cm' 2 s~ 1 , corresponding to an afterglow luminosity of 
2.8 x 10 43 erg sT l using the redshift (z = 0.53) determined from 
the thermal plasma model fit to the data (see Sect. 3.1). The 
second source is fainter (1.2 x 10~ 14 ergcrrr 2 s~') and shows 
no evidence for variability. 

2.2. GRB010220 

GRB 010220 was detected by BeppoSAX at 22:51:07 UT on 
20 February 2001 and subsequently localised to within a 4' ra- 
dius circle; the burst duration was ~ 40 s, with ~ 660 count s _1 
peak flux in the 40-700keV band (Manzo et al. 2001). No 
optical counterpart was detected to a limiting magnitude of 
R^ 23.5 (Berger et al. 2001), however the burst position is 
in the Galactic plane (b n = 1.41), where the Galactic nebula 
IC 1805 lies along the line of sight (Castro-Tirado et al. 2001). 
The Galactic hydrogen absorbing column is 8.6 x 10 21 crrT 2 
(Dickey & Lockman 1990, using the FTOOL nh). At 14.8 
hours after the burst, XMM-Newton began observing at the 
coordinates of the BeppoSAX error-circle. The observation of 
GRB 010220 was contaminated by a high and variable back- 
ground. To mitigate this effect, only data where the background 
rate was relatively low were used. The EPIC-pn exposure was 
43 ks; screening for flares left 20 ks of good data. Full frame 
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mode and the medium filter were used for the observation. No 
data were available from the EPIC-MOS cameras. 

The data were processed and reduced with the SAS ver- 
sion 5.2. A source extraction region of 20" radius and an off- 
source background extraction region of 80" radius were cho- 
sen. Both single and double pattern, good (FLAG=0) events 
were used with the ready-made response matrix 1 provided by 
the XMM-Newton SOC. The spectra were binned with a min- 
imum of twenty counts per bin. Source positions were deter- 
mined as noted above. 

Four sources were detected in the BeppoSAX error circle; 
the brightest source is 18" from the centre, with the others at 
> 3'. No variability is detected in any of the sources at greater 
than 2cr confidence. We assume that the brightest source, at 
coordinates (J2000): R.A. 02h37m0 1.66s, Dec. +61°45'56.01" 
with a 68% confidence error circle of 1.2" radius is the GRB 
afterglow. This source decays at 94% confidence, with a best-fit 
decay index,/? = 1.2 + 1.0. Its mean 0.2-10.0keV flux was 3.3x 
10" 14 ergcm" 2 s _1 (corresponding to an afterglow luminosity 
of 7.3 x 10 43 erg s" 1 in the rest frame 0.2-10.0keV band based 
on the best-fit redshift (z = 1 .0) from the thermal plasma fit to 
the data). 

3. Spectral fitting 

In order to test the wider applicability of the collisionally- 
ionised plasma model proposed by Reeves et al. (2002) to 
explain the XMM-Newton observations of the afterglow of 
GRB 011211, our new GRB X-ray afterglow data were fit with 
the same set of models. They were a) a power-law, b) a power- 
law with a variable, cold, redshifted absorber, c) a variable 
abundance collisionally-ionised plasma model (the MEKAL 
model, Mewe et al. 1985, Liedahl et al. 1995) with the abun- 
dances of Mg, Si, S, Ar and Ca fit jointly, Ni allowed to vary 
freely and all other elements fixed at the solar value and d) 
an ionised reflection model (Ballantyne & Ramirez-Ruiz 2001) 
where the emission arises purely from the X-ray flux scattered 
off material with twice the solar elemental abundance, also 
modified by a cold, redshifted absorber. In all cases an absorber 
fixed at the Galactic value was included. 

In order to test the significance of the fit improvements, 
10000 spectra were simulated using the parameters derived 
from fitting an absorbed power-law (model b) to the data in 
each case. Models b) and c) were fit to the simulated data and 
the difference in^ 2 computed. The results of these tests are re- 
ported in Table 1 . Using a random gaussian distribution of ini- 
tial parameters centred on the parameters derived from fitting 
an absorbed power-law to the data, a second set of 10 000 sim- 
ulated spectra was generated. This second set of simulations 
yielded results consistent with the first. 

3.1. GRB 001 025 A 

The spectrum for this afterglow is significantly better fit 
by the thermal plasma model than the absorbed power-law 
(Table 1) — only 13 in 10000 simulations are as good a fit. 
While the continuum shape dictates that the plasma tempera- 
tures in the simulated spectra will be consistent with each other 
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Table 1. ^/degrees of freedom for model fits to both GRB afterglows. 
Model a) is a power-law. Model b), is the same as model a) modified 
by a cold, redshifted absorber, c) is a variable abundance collision- 
ally ionised plasma model with the abundances of Mg, Si, S, Ar and 
Ca fit jointly, Ni allowed to vary freely and all other elements fixed 
at the solar value. Model d) is an ionised reflection model where the 
emission arises purely from the X-ray flux scattered off material with 
twice the solar elemental abundance, also modified by a cold, red- 
shifted absorber. In all cases an absorber fixed at the Galactic value 
was included. Column four is the percentage of spectral fits (simu- 
lated using model b) with greater x 1 differences between models b) 
and c) than the real data. 



GRB 


a 


b 


MC Prob. 


c 


d 


00 1025 A 
010220 


95.7/66 
15.4/15 


68.6/64 
14.8/14 


0.13% 
0.16% 


54.1/63 
5.2/13 


65.5/63 
6.7/12 



Table 2. Best-fit parameters from spectral fits to the afterglow data. 
In the upper section, parameters of the fit to the collisionally ionised 
plasma model (model c in Sect. 3) are presented. Parameters of the 
best fit to an absorbed power-law (model b) are included for reference 
in the lower section. Values in parentheses are 90% confidence limits 
for one interesting parameter. 



GRB 


T (keV) 


z 


Mg/Si/S/Ar/Ca 


Ni 








(Solar Abundances) 


001025A 


3.4 


0.53 


2.0 


26 




(2.9-3.9) 


(0.50-0.55) 


(0.7-3.7) 


(14-40) 


010220 


6.0 


1.0 


10 


95 




(3.6-12.7) 


(0.97-1.07) 


frozen 


(32-248) 




r 


z 


N H (10 21 cm" 2 ) 




001025A 


2.50 


0.5 


3 






(2.26-2.97) 


(0-7.1) 


(1-184) 




010220 


2.1 


1.0 


16 






(1.5-3.1) 


frozen 


(0-67) 





and with the fit to the real data, only four of the redshifts de- 
termined from the thirteen best-fit simulated spectra are within 
90% confidence limits of the value from the real data, reinforc- 
ing the fact that the emission lines are not systematic devia- 
tions. 

The need for the thermal emission component arises from 
the soft excess observed between ~ 0.5 and 2keV (Fig. 1), 
which we suggest may be due to the blend of lines from Mg, 
Si, S, Ar and Ni-L. The thermal plasma model fit to the data 
yields a temperature of 3.4 keV and a redshift of 0.53 and does 
not require super-solar abundance of light metals (Table 2). 

In order to parameterise the line emission, a power-law 
with Galactic absorption and Gaussian emission lines was fit 
to the data, allowing the line widths to vary together. For five 
emission lines, the x 1 is 40. 1 for 55 degrees of freedom, (giv- 
ing a null hypothesis probability of 3 x 10~ 4 (3.6 cr) over the 
absorbed power-law and 5 x 10~ 7 (5.0 cr) over a power-law 
with Galactic absorption. Adding lines to a bremsstrahlung 
continuum yielded similar significances, (3 x 10~ 4 (3.6 cr) 
over the absorbed bremsstrahlung and 2 x 10~ 4 (3.8 cr) over 
a bremsstrahlung with Galactic absorption), though the fit was 
worse, x 2 = 45.4 for 55 degrees of freedom. 



The best-fit energies for the lines are: 4.7^4, 2.2 ±0.1, 
1.64 ± 0.07, 1.16 + 0.05 and 0.80^ keV with a FWHM of 
240^0 eV with individual line significances (null hypothesis 
probabilities) of 0.12, 0.02, 8 x 10~ 4 , 8 x 10~ 4 and 4 x 10~ 4 
respectively. Single line significances were determined by us- 
ing the best-fit Galactic-absorbed power-law with five Gaussian 
emission lines and removing individual lines to assess their f- 
statistic. Probabilities were then determined using the f-test. 
The observed equivalent width is ~ 800 eV for the 4.7 keV line 
and ~ 400 eV for the others. Fitting the Ka emission lines of 
Mgxn (1.46keV), Sixiv (1.99keV), Sxvi (2.60 keV), Arxvm 
(3.30keV) and 56 Nixxvm (8.10keV), gives a best-fit redshift 
of 0.7+g j. This redshift is within two standard deviations of 
the redshift determined from the thermal plasma fit (Table 2). 
Allowing for an outflow velocity of 0.0 - 0.1 c (Reeves et al. 
2002), we conclude that the redshift of the host galaxy is likely 
to lie in the range 0.5-1.2. 

It has been suggested that reflection of synchrotron emis- 
sion from ionised material could produce the emission lines 
observed in GRB afterglows (e.g. Ballantyne et al. 2002). An 
ionised reflection model is as good a fit to the data as an ab- 
sorbed power-law; however, adding a power-law continuum to 
the reflection component disimproves the fit significantly, in- 
dicating that the afterglow must be viewed so that little of the 
continuum is observed, if this fit is to be acceptable. This im- 
plies a line-of-sight outside the GRB cone, a possibility ruled 
out by the detection of the GRB in the first place. In any case, 
the ionised pure reflection model does not fit the soft excess as 
well as the thermal plasma model (Table 1). Modifications to 
the reflection model (such as including contributions from S, 
Ar and Ca) could improve the fit to the data. 

3.2. GRB 01 0220 

A power-law model with Galactic absorption is an acceptable 
fit to the data, however the single clearest deviation in the spec- 
trum of GRB 010220 (Fig. 2) is the feature near 3.9 keV; adding 
an unresolved Gaussian line (equivalent width = 1.8^° ^keV) 
to the power-law fit improves the fit at > 99% significance 
(X 2 = 7.5 for 13 degrees of freedom) and is the only emission 
feature in the spectrum. While the relatively poor statistics do 
not enable us to detect individual emission lines at lower ener- 
gies (e.g. Si or S Ka), the thermal plasma model is significantly 
preferred to the absorbed power-law model with a null hypoth- 
esis probability of 0.0016 using statistics from the Monte Carlo 
simulations. Allowing the abundance of Fe to vary did not sig- 
nificantly improve the thermal plasma fit as the energy of the 
Fe Ka emission is too low relative to the emission at low ener- 
gies and is not broad enough to fit the excess alone. The best-fit 
Fe abundance is consistent with solar abundance; allowing Ni 
to vary improves the fit at > 99% confidence as it fits most of 
the 3.9 keV emission feature. 

The pure ionised reflection model in this case is a better fit 
than the absorbed power-law model (Table 1) primarily due to 
fitting Fe emission (at lower redshift than the thermal plasma 
model) to the ~ 3.9 keV feature, however the addition of any 
power-law continuum disimproves the fit as for GRB 00 1025 A. 
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Fig. 1. EPIC-pn (crosses) and combined MOS (filled triangles) spec- 
trum of the afterglow of GRB001025A; a) data fit to a Galactic- 
absorbed power-law model; b) residuals to plot a); c) fit residuals to a 
Galactic-absorbed power-law model with variable redshifted absorp- 
tion; d) fit residuals to a collisionally-ionised plasma model with vari- 
able abundance of Mg, Si, S, Ar and Ca and freely variable Ni abun- 
dance (VMEKAL), absorbed by the Galactic column; e) fit residuals 
to a Galactic-absorbed power-law and five Gaussian emission lines 
model, allowing the line widths to vary jointly. 



4. Discussion 

It is clear that for the X-ray afterglows of GRBs 001025A and 
010220, the type of collisionally-ionised plasma model fit to 
the afterglow of GRB 011211 (Reeves et al. 2002) is a better 
fit to the data than an absorbed power-law. The detection of a 
thermal spectrum in three GRB afterglows with XMM-Newton 
indicates that this behaviour is relatively common hours to days 
after the burst among the long duration GRBs. Thermal emis- 
sion may therefore make a significant contribution to the total 
afterglow luminosity of a GRB. It is worth noting that these 
are among the best quality soft X-ray spectra of GRB after- 
glows recorded to date and the absence of detection of such 




chonnel energy (keV) 



Fig. 2. EPIC-pn spectrum of the afterglow of GRB 010220, 20 ks ex- 
posure, fit with a power-law with Galactic absorption. 



thermal spectra with BeppoSAX is not suprising given XMM- 
Newton^ much greater sensitivity and better spectral resolu- 
tion and the relatively high low-energy cut-off of the MECS in- 
strument (~ 2 keV). A re-examination of the ASCA, BeppoSAX 
and Chandra data of the afterglows of GRBs 970508, 970828 
(where the 5.04 keV feature could be caused by blueshifted 
56 Ni Ka emission blueshifted by 0.2c), 991216, 990705 and 
000214 may be worthwhile in the light of these results. 

In none of the three cases (GRBs 00 1025 A, 010220 and 
011211) is an excess abundance of Fe required and in both af- 
terglow spectra reported here a large over-abundance of Ni (or 
possibly Co) is required, consistent with the idea that the SN- 
GRB lag must be of the order of days, not months, though it 
is worth noting that this timescale can be extended where the 
material is ionised (McLaughlin & Wijers 2002). 

Using the model proposed by Reeves et al. (2002) to ex- 
plain the observations of GRB 01 1211 it is difficult to place 
precise limits on the size of the SN shell, the expansion speed 
and the time delay between the burst and the initial SN from the 
spectra of these afterglows as the redshifts of the host galaxies 
are unknown and the redshifts derived from the X-ray spectra 
are not as tightly constrained as in GRB 01 1211. It is possi- 
ble however, to constrain the minimum radius of the SN ejecta 
shell and its mass on the assumption that the duration of the 
thermal emission is due primarily to light-crossing time over 
the illuminated shell. With no collimation of the GRB and 
where the thermal emission lasts at least as long as the time 
between the burst and the afterglow observation (~ 7 and 30 
hours), the minimum shell radii are 8 x 10 14 cm and 3 x 10 15 cm 
corresponding to SN-GRB delays of at least 3 and 12 days for 
GRB 010220 and GRB 00 1025 A respectively and illuminated 
shell masses of 0.9 and 1 . 1 Mq using the best fit emission 
measure, temperature and redshifts from the MEKAL model. 
Reeves et al. (2002) calculate a delay of a few days and at 
most two weeks between the SN and the GRB, consistent with 
these results, though the SN-GRB time delay distribution has 
so far not been defined. The time delay for GRB 00 1025 A is 
suggestive since the half-life of the most abundant product of 
Si burning, 56 Ni, is 6.1 days decaying primarily to 56 Co with 
a half-life of 77.3 days to 56 Fe (although see McLaughlin & 
Wijers 2002, on half-lives of ionised species). In this scenario, 
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much of the emission observed near 1 keV in the spectrum 
of GRB 00 1025 A could then be due to L-shell emission from 
56 Co. The spectrum however is not sufficiently good to investi- 
gate this possibility. 

Recent analysis of the XMM-Newton spectrum of 
GRB 020322, with a signal to noise ratio comparable to the 
X-ray spectrum of GRB 01 1211, shows no evidence of emis- 
sion features (Watson et al. in prep.) 15 hours after the burst; 
evidence that if thermal emission occurs in every afterglow, it 
may be relatively short-lived. 

The prospect that thermal emission may contribute signifi- 
cantly to the X-ray luminosity of afterglows could resolve the 
problem of bursts with apparently high X-ray column densities 
but low optical extinctions, 10-100 times smaller than expected 
(Galama & Wijers 2001). The curvature of low-quality thermal 
X-ray spectra can be approximated by an absorbed power-law, 
giving much larger absorbing columns than are required by the 
thermal spectra. 

Another implication of the data reported here is that GRBs 
without significant optical emission, "dark GRBs" (Reichart & 
Yost 2001, Djorgovski et al. 2001), may not always need to 
have significant extinction (Taylor et al. 2000, Reichart & Yost 
2001, Piro et al. 2002), be intrinsically different (Lazzati et al. 
2002) or be very distant (Lamb & Reichart 2000). A soft X-ray 
spectrum dominated by thermal emission will lead to an over- 
estimation of the expected optical flux where this is predicted 
assuming a synchrotron-dominated X-ray flux below lOkeV 
(Djorgovski etal. 2001). 

5. Conclusions 

XMM-Newton detected the afterglows of GRB 00 1025 A and 
GRB 010220 the former showing a decaying lightcurve. 
Positions accurate to ~ 1" were determined; this makes it fea- 
sible for optical spectroscopy to be used to confirm the X-ray 
redshift reported in this paper for the host galaxies of these 
GRBs, in particular for the host of GRB 001025 A, where the 
extinction is relatively low. In both cases, the X-ray spectra 
are significantly better fit by the thermal plasma model pro- 
posed by Reeves et al. (2002) to explain the line features in 
GRB 01 1211 than by an absorbed power-law and in the case 
of GRB 00 1025 A the thermal plasma model is a significantly 
better fit than an ionised reflection model. The parameters de- 
termined from these thermal fits are consistent with their sug- 
gested scenario. It seems likely that a large contribution from 
highly-ionised light metals is a common feature in the X-ray 
spectra of GRB afterglows hours to days after the burst and 
the detection of thermal emission in three of the four XMM- 
New ton-detected afterglows implies that this may be a signif- 
icant component in the total afterglow luminosity of all long 
GRBs. 
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